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The value of the V^d matrix element of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix can be derived from nuclear superallowed beta decays, neutron decay, and 
pion beta decay. We survey current world data for all three. Today, the most 
precise value of y^^j comes from the nuclear decays; however, the precision is 
limited not by experimental error but by the estimated uncertainty in theoretical 
corrections. Experimental uncertainty does limit the neutron-decay result, which, 
though statistically consistent with the nuclear result, is approximately a factor of 
three poorer in precision. The value obtained for V„d leads to a result that differs 
at the 98% confidence level from the unitarity condition for the CKM matrix. We 
examine the reliability of the small calculated corrections that have been applied 
to the data, and assess the likelihood of even higher quality nuclear data becoming 
available to confirm or deny the discrepancy. Some of the required experiments 
depend upon the availability of intense radioactive beams. Others are possible 
today. 



1 Introduction 

The Cabibbo-Kobayashi-Maskawa matrix'^'^ relates the quark eigenstates of 
the weak interaction with the quark mass eigenstates (unprimed) 

'Vud Vus Vub\ fd\ 
Vcd Vcs Kb s (1) 
Vtd Vts Vtb) \bj 

and, as such, the matrix is unitary. Thus there are many relationships among 
the nine elements of the matrix that can be tested by experiment. The leading 
element, Vud, only depends on quarks in the first generation and so is the 
element that can be determined most precisely. Here we will discuss the current 
status of Vud and the unitarity test as it relates to the elements in the first 
row: 



Vud ~^ Vus ^~ Vub ^' 



(2) 



In examining the unitarity test, we will adopt the Particle Data Group 
(PDG02) ^ recommendations for Vus and Vub- Wc note that for Vus PDG02 
recommends only the value determined from K^s decay, | Vus |= 0.2196±0.0023, 
arguing that the value obtained from hyperon decays suffers from theoretical 
uncertainties due to first-order SU(3) symmetry-breaking effects in the axial- 
vector couplings. It now appears, though, from preliminary results that the 
experimental result for the e3 branch of the decay may also be in doubt. 
For the time being, we continue to use the value of Vus from PDG02 but 
acknowledge that this result may have to change in future. As to Vub, its value 
is small, | Vub |= 0.0036 ± 0.0010, and consequently it has a negligible impact 
on the unitarity test, Eq. (2). 

2 The value of Vud 

The value of Vud can be determined from three distinct sources: nuclear super- 
allowed Fermi beta decays, the decay of the free neutron, and pion beta decay. 
We discuss each in turn. 

2.1 Nuclear superallowed Fermi beta decays 

Nuclei have the singular advantage that transitions with spcicific characteristics 
can be selected and then isolated for study. One example is the superallowed 
0+ 0+ beta transitions, which depend uniquely on the vector part of the 
weak interaction. Furthermore, in the allowed approximation, the nuclear 
matrix element for these transitions is given by the expectation value of the 
isospin ladder operator, which is independent of any details of nuclear struc- 
ture and is given simply as an SU(2) Clebsch-Gordan coefficient. Thus, the 
experimentally determined /^-values are expected to be very nearly the same 
for all 0+ — > 0+ transitions between states of a particular isospin, regardless 
of the nuclei involved. Naturally, there are corrections to this simple picture 
coming from electromagnetic effects, but these corrections are small - of or- 
der 1% - and calculable. Thus, if we write Sr as the nucleus-dependent part 
of the radiative correction, as the nucleus-independent part of the radia- 
tive correction, and 5c as the isospin symmetry-breaking correction, then the 
experimental ft-value can be expressed as follows: 

ft{l + 6R){l-Sc)=J't= (f+ = constant, (3) 

wherein = 271^ \n 2h{hcf/{meC^)^, which has the value K/{hc)^ = (8120.271± 
0.012) X 10~^° GeV~^s; and Gy is the weak vector coupling constant: Gy = 
GfVud, with Gp being the weak coupling constant from purely leptonic muon 



Table 1: Experimental results {Qec, ti/2 ^■iid branching ratio, R) and calculated correction, 
PeCi foi^ 0+ — » 0+ transitions. The other calculated corrections, 5ji and 5c, are given in 
Tables 3 and 4 respectively. 





Qec 


*l/2 


R. 


Pec 


ft 


rt 




(keV) 


(ms) 


(%) 


(%) 


(s) 


(s) 


IOC 


1907.77(9) 


19290(12) 


1.4645(19) 


0.296 


3038.7(45) 


3072.9(48) 


14 o 


2830.51(22) 


70603(18) 


99.336(10) 


0.087 


,3038.1(18) 


,3069.4(26) 




4232.42(35) 


6.344.9(19) 


> 99.97 


0.083 


.3035.8(17) 


3071.4(22) 


34 Cl 


5491.71(22) 


1525.76(88) 


> 99.988 


0.078 


.3048.4(19) 


3070.6(25) 


38m 


6044.34(12) 


923.95(64) 


> 99.998 


0.082 


.3049.5(21) 


3070.9(27) 


42 Sc 


6425.58(28) 


680.72(26) 


99.9941(14) 


0.095 


3045.1(14) 


3075.7(24) 


46v 


7050.63(69) 


422.51(11) 


99.9848(13) 


0.096 


3044.6(18) 


,3074.4(27) 


50 Mn 


7632.39(28) 


283.25(14) 


99.942(3) 


0.100 


3043.7(16) 


,3072.9(28) 


54 Co 


8242.56(28) 


193.270(63) 


99.9955(6) 


0.104 


3045.8(11) 


3072.1(27) 










Average, Tt 


3072.2(8) 
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decay. Thus, to extract Vud from experimental data, tlic^ procedure is to de- 
termine the JTf-vahies for a variety of different nuclei having the same isospin, 
and then to test if they are self-consistent. If they are, their average is used to 
determine a value for Gv and, from it, Vud- 

It is now convenient to separate the radiative correction into two terms: 

5R = 5'ii + 5ns (4) 

where the first term, 5^, is a function of the electron's maximum energy and 
the charge of the daughter nucleus, Z; it therefore depends on the particular 
nuclear decay, but is independent of nuclear structure. To emphasize the dif- 
ferent sensitivities of the correction terms, we re- write the expression for Tt 
as 

J't = ft{l + 5'j,){l + 5NS-5c) (5) 

where the first correction term in brackets is independent of nuclear structure, 
while the second incorporates the structure-dependent terms. 

To date, superallowed 0+ 0"*" transitions have been measured to ±0.1% 
precision or better in the decays of nine nuclei ranging from to ^'*Co. 
World data on Q-values, lifetimes and branching ratios - the results of over 
100 independent measurements - were thoroughly surveyed^ in 1989 and then 
updated several times since then, most recently at the time of the last WEIN 
conference^. These results appear in the first three columns of Table 1. Using 
the calculated electron-capture probabilities ^ given in the next column, we 
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Figure 1: ^t-valucs for the nine precision data, and the best least-squares one-parameter fit. 



obtain the "uncorrected" /t-values listed in column 5 with partial half-lives 
determined from the formula t = ^1/2(1 + Pec)/R- 

We save a detailed discussion of the radiative and isospin symmetry- 
breaking corrections for Sec. 3. In the present context, though, it should be 
noted that the values for and 5c used to derive J-t were taken from the last 
column of Tables 3 and 4, respectively; each value results from more than one 
independent calculation. In the case of 6r, the calculations are in complete 
accord with one another; for 5c ^ we have used an average of two indepen- 
dent calculations with assigned uncertainties that reflect the (small) scatter 
between them. Thus, in a real sense, both experimentally and theoretically, 
the .Ft-values given in Table 1 and plotted in Fig. 1 represent the totality of 
current world knowledge. The uncertainties reflect the experimental uncer- 
tainties and an estimate of the relative theoretical uncertainties in 5c There is 
no statistically significant evidence of inconsistencies in the data jv = 1.1), 
thus verifying the expectation of CVC at the level of 3 x 10^*, the fractional 
uncertainty quoted on the average ^i-value. 

In using this average J^t- value to determine Gv, we must account for ad- 
ditional uncertainty: viz 



:'^t = 3072.3 ±0.9 ±1.1, (6) 

where the first error is the statistical error of the fit (as illustrated in Fig. 1), 
and the second is an error related to the systematic difference between the two 
calculations of Sc by Towner and Hardy ^ and by Ormand and Brown ^ that 
we have combined in reaching this result. (For a more complete discussion of 
how we treat these theoretical uncertainties, see reference^.) We now add the 
two errors linearly to obtain the value we use in subsequent analysis: 

^= 3072.3 ±2.0. (7) 

The value of Vud is obtained by relating the vector constant, Gv, deter- 
mined from this J^t value, to the weak coupling constant from muon decay, 
G^/{ncf = (1.16639 ± 0.00001) x 10"^ GeV-^, according to: 

T/-2 ^ ^ /ON 

"■^ 2G2(l + Ax)jri' ^ ^ 

The result obtained is 

I Vud 1= 0.9740 ± 0.0005, [Nuclear] (9) 

where the nucleus-independent radiative correction has been set at 

= (2.40 ± 0.08)%. (10) 

Note this value differs slightly (but within errors) from an earlier value 
because of the decision by Sirlin-'^^ to centre the cut-off parameter TOa, where 
(mai/2) < TOa < 2mai, exactly at the ai-meson mass when evaluating the 
axial contribution to the radiative-correction loop graph. 

Prom the value of Vud given in Eq. (9), the unitarity sum, Eq. (2), becomes 



0.9968 ± 0.0014, [Nuclear] (11) 



which fails to meet unity by 2.2 standard deviations. In connection with this 
result, we note the following two points: 

(a) The error bar associated with | Vud \ in Eq. (9) is not predominantly 
experimental in origin. In fact, if experiment were the sole contributor, the 
uncertainty would be only ±0.0001. The largest contributions to the | Vud \ 
error bar come from (±0.0004) and Sc (±0.0003). 

(b) The unitarity result in Eq. (11) depends on the values of nuclear 
structure-dependent corrections. In Sec. 3 we will examine whether the failure 



to meet unitarity can be repaired by reasonable adjustments to these correc- 
tions. Our conclusion is largely negative. Other speculative possibilities are 
presented in Sec. 4. 

2.2 Neutron decay 

On the one hand, free neutron decay has an advantage over nuclear decays 
sinc;c thcirc are no nuclear-structure dependent corrections to be calculated. 
On the other hand, it has the disadvantage that it is not purely vector-like 
but has a mix of vector and axial-vector contributions. Thus, in addition to 
a lifetime measurement, a correlation experiment is also required to separate 
the vector and axial-vector pieces. Both types of experiment present serious 
experimental challenges. The value of Vud is determined from the expression 



where A is the ratio of axial-vector and vector effective coupling constants. 



A = G'JG'-^, with Gf = (1 + A*) and G'-^ = (1 + A^). Here A* and 



arc the nucleus-independent radiative corrections. With the experimental 
measurement of A being a determination of G'^/G'^, the actual value of A* is 
not required for the evaluation of V^^. In Eq. (12), / is the statistical rate 
function and (5/?, the nucleus-dependent radiative correction evaluated for the 
case of a neutron. Wilkinson has evaluated the product f{l + 6R). His value 
was revised by Towner and Hardy who incorporated the current best Q- value 
to obtain /(I + 6r) = 1.71489 ± 0.00002. Lastly, t„ is the mean lifetime for 
neutron decay. 

A survey of world data on neutron decay appears in Table 2. The life- 
time measurements, when considered as a single body of data, are statistically 
consistent. However the measurements of A are not jv = 2.6) and, as a 
consequence, the uncertainty quoted in the table for the overall average value 
of A has been scaled by a factor of 1.6. Inserting the average values for t„ and 
A into Eq. (12), we determine the value of Vud to be 




Kj In 2 



(12) 



G|(l + A-)(l + 3A2)/(l + 5;j)r„' 



\Vud\= 0.9745 ± 0.0016, 



[Neutron] 



(13) 



and the unitarity sum to be 



^ y„2 ^ 9978 _^ 0.0033, 

i 



[Neutron] 



(14) 



Table 2: Experimental results" for neutron decay 



Ak-Uicjd 


Measure 






T„(s) n beam 


891 ±9" 






n trap 


877 ± 10 15 


887.6 ± 3.0 1^ 






888.4 ±3.31^ 


882.6 ± 2.7 1* 






889.2 ± 4.8 1^ 


885.4 ± 1.0 20 




Overall Average 






885.6 ± 0.8 


A P-asym. 


-1.254 ±0.015^ 21 


-1.257 ±0.012'' 22 






-1.262 ±0.005 23 


-1.2594 ± 0.0038 2" 






-1.266 ±0.00425 


-1.2739 ± 0.0019 26 




(3 and v asym 


-1.2686 ±0.0047 27 






corr. 


-1.259 ±0.01728 






Overall Average 






-1.2690 ± 0.0022 



" Following the practice used for superallowed decay, we retain only those measurements 
with uncertainties that are within a factor of ten of the most precise measurement for each 
quantity. All such measurements, of which we are aware, that have not been withdrawn by 
their authors are listed. 



Corrected for weak magnetism and recoil following ref. 



a value that agrees with unitarity and with the nuclear result, Eq. (11) which 
is over a factor of two more precise. In connection with this result, we note 
the following two points: 

(a) For neutron decay, the error bar associated with | Vud \ in Eq- (13) is 
some three times larger than the error bar obtained from nuclear decays, Eq. 
(9); however, in contrast with the latter case, it is predominantly experimental 
in origin. The largest theoretical contribution to the | Vud \ error bar comes (at 
the level of ±0.0004) from A^, a correction that is common to both nuclear 
and neutron decays. 

(b) Currently, the theoretical uncertainty on dominates the nuclear 
result for | Vud \- As experimental results for the neutron improve, AJJ will 
eventually dominate the neutron result too. Therefore, so long as AJ( remains 
at its current level of uncertainty, the neutron results will never be able to test 
unitarity with significantly better precision than the nuclear decays do now, 
in spite of their independence of 6c- They will, of course, be able to provide 
an important test of whether there are some systematic problems with the 
nuclear-dependent corrections, which are not now anticipated in the theoretical 
uncertainty quoted in Eq. (9). 



2.3 Pion beta decay 

Like neutron decay, pion beta decay has an advantage over nuclear decays 
in that there are no nuclear structure-dependent corrections to be made. It 
also has the same advantage as the nuclear decays in being a purely vector 
transition, in its case 0~ — > 0~, so no separation of vector and axial- vector 
components is required. Its major disadvantage, however, is that pion beta 
decay, 7r+ 7r°e+fe, is a very weak branch, of the order of 10~^. This results 
in severe experimental limitations. For the pion beta decay, the value of Vud 
is determined from the expression 

y2 ^ K/\n2 

"-^ GUl + Al)M2f{l + 6n)rm/BR' ^ 

where / is the approximate statistical rate function 

1 / A ^ ^ 



with A being the pion mass difference (A = m„+ — m^o ) and , the electron 
mass. The mass difference is known with high precision from the work of 
Crawford et al?^ to be A = (4.5936 ± 0.0005) MeV. The factors /i and /2 are 
corrections to /, and are easily calculable hmctions""^ of A/m.j^+ with values of 
/i = 0.941039 and /2 = 0.951439. The nucleus-dependent radiative correction 
evaluated for the case of the pion is Sn = (1.05 ± 0.15)%. Finally, Tm and BR 
are the pion mean lifetime and branching ratio, respectively. 

Two precise lifetime measurements '^^''^^ were published in 1995 and the 
PDG02 average is 

Tm = (2.6033 ± 0.0005) x 10~^ s. (17) 
The branching ratio is principally from McFarlane et al?'^: 

BR = (1.025 ± 0.034) x 10"^. (18) 
Inserting Eqs. (17) and (18) into Eq. (15), we determine the value of Vud to be 

1 14^1= 0.9670 ±0.0161, [Pion] (19) 

and the unitarity sum 



^y^^i =0.9833 ±0.0311, 

i 



[Pion] 



(20) 



satisfying the unitarity condition but with comparativefy large uncertainty. 

The error on | Vud \ is entirely due to the uncertainty in the pion branching 
ratio. There is an experiment currently underway at PSI designed to improve 
this branching ratio by a factor of eight. If successful, the error on | Vud | would 
be reduced to ±0.0022. still less precise than the current results from both 
neutron and nuclear decays. Of course, ultimately it too will be limited by the 
theoretical uncertainty on A^. 

3 Theoretical corrections in nuclear decays 

In Sec. 2.1 we noted that the value of Vud determined from nuclear decays 
- the most precise result available - resulted in the unitarity test among the 
elements of the first row of the CKM matrix not being satisfied by two stan- 
dard deviations. Here, we discuss the theoretical corrections involved in the 
determination, and assess whether the failure to meet unitarity can be re- 
moved by reasonable adjustments in these calculations. To restore unitarity, 
the calculated radiative corrections for all nine nuclear transitions would have 
to be shifted downwards by 0.3% {i.e. as much as one-quarter of their current 
value), or the calculated isospin symmetry-breaking correction shifted upwards 
by 0.3% (over one- half their value), or some combination of the two. 

3.1 Radiative corrections 

As mentioned in Sec. 2.1, the radiative correction is conveniently divided into 
terms that are nucleus-dependent, Sn, and terms that are not, A^. These are 
written 



Sr = 6'ij^ + Sns 

S'n = ^[9iEm) + S2 + Ss] 

cx 

= 7^ [4 ln(mz/TOp) + ln(mp/mA) + 2CBor„] + • • • , (21) 

ZTT 

where the ellipses represent further small terms of order 0.1%. In these equa- 
tions, is the maximum electron energy in beta decay, ruz the Z-boson 
mass, rriA the ai-nicson mass, and S2 and 6^ the oidci-Za^ and -Z^a^ contri- 
butions respectively. The function g{E,,, Em), which is a function of electron 
energy, was first defined by Sirlin '^'^ as part of the order-a universal photonic 
contribution arising from the; weak vector current; it is here averaged over 
the electron spectrum to give 'g{Em). Finally, the terms Ceorn and Sns come 
from the order-a axial-vector photonic contributions: the former accounts for 



Table 3; Calculated nucleus-dependent radiative correction, Sn, in percent units, and the 
component contributions as identified in Eq. (21). 















10c 


1.468 


0.182 


0.005 


1.65(1) 


-0.360(35) 


140 


1.286 


0.227 


0.008 


1.52(1) 


-0.250(50) 


26m ^2 


1.110 


0.325 


0.021 


1.46(2) 


0.009(20) 


34Q 


1.002 


0.388 


0.034 


1.42(3) 


-0.085(15) 


38m 


0.964 


0.413 


0.041 


1.42(4) 


-0.100(15) 




0.939 


0.448 


0.049 


1.44(5) 


0.030(20) 


46v 


0.903 


0.468 


0.057 


1.43(6) 


-0.040(7) 


50Mn 


0.873 


0.494 


0.065 


1.43(7) 


-0.042(7) 


54C0 


0.843 


0.507 


0.073 


1.42(7) 


-0.029(7) 



single- nuclcon contributions, while the latter covers two-nucleon contributions 
and is consequently dependent on nuclear structure. 

Calculated values for all four components of 6n are given in Table 3. There 
have been two independent calculations of both S2 and ^3; they arc com- 
pletely consistent with one another if proper account is taken of finite-size 
effects in the nuclear charge distribution. The values listed in Table 3 are 
our recalculations ^ using the formulas of Sirlin "^^ but incorporating a Fermi 
charge-density distribution for the nucleus. Note that we have followed Sirlin 
in assigning an uncertainty equal to {a/27T)S3 as an estimate of the error made 
in stopping the calculation at that order. Also appearing in the table arc our 
adopted values for Sns, which we have taken from our recent re-evaluation^. 

In assessing the changes in Sr that would be required in order to restore 
unitarity, it is helpful to rewrite Eq. (21) in terms of the typical values taken 
by its components: viz 

5r ~ 1.00 + 0.40 + 0.05. (22) 

Thus, ~ 1.4%. If the failure to obtain unitarity in the CKM matrix with 
Vud from nuclear beta decay is due to the value of J^, then all S'j^ values must 
be reduced to ~ 1.1%. This is not likely. The leading term, 1.00%, involves 
standard QED and is well verified. The order-Za^ term, 0.40%, while less 
secure has been calculated twice ^''■^^ independently, with results in accord. 

Taking a similar approach for the nucleus- independent radiative correction, 
we write 



A^ = 2.12 -0.03 + 0.20 + 0.1% 2.4%, 



(23) 



of which the first term, the leading logarithm, is imambiguous. Again, to 
achieve unitarity of the CKM matrix, would have to be reduced to 2.1%: 
i.e. all terms other than the leading logarithm must sum to zero. This also 
seems unlikely. 

3.2 Isospin symmetry-breaking corrections 

Because the leading terms in the radiative corrections are so well founded, 
attention has focussed more on the isospin symmetry-breaking correction. Al- 
though smaller than the radiative correction, the isospin symmetry-breaking 
correction is clearly sensitive to nuclear-structure issues. It comes about be- 
cause Coulomb and charge-dependent nuclear forces destroy isospin symmetry 
between the initial and final states in superallowed beta-decay. The conse- 
quences are twofold: there are different degrees of configuration mixing in 
the two states, and, because their binding energies are not identical, their 
radial wave functions differ. Thus, we accommodate^ both effects by writing 
= Sci+Sc2- Constraints can be placed on the calcidation of 6ci by insisting 
that it reproduce the measured coefficients of the isobaric mass multiplet equa- 
tion. Constraints arc placed on Sc2 by insisting that the asymptotic forms of 
the proton and neutron radial functions match the known separation energies. 

The results of several calculations for 5c are shown in Table 4. The values 
in the first column are those calculated by the methods developed by Towner, 
Hardy and Harvey^ and refined in more recent publications They result 
from shell-model calculations to determine Sci, and full-parentage expansions 
in terms of Woods-Saxon radial wave functions to obtain 6c2- Ormand and 
Brown, whose values^ for Sc appear in column 2, also employed the shell model 
for calculating 5ci but derived Sc2 from a self-consistent Hartree-Fock calcu- 
lation. Both of these independent calculations for Sc reproduce the measured 
coefficients of the relevant isobaric multiplet mass equation, the known pro- 
ton and neutron separation energies, and the measured /t-values of the weak 
non-analogue 0+ — > O"'" transitions'*" where they are known. In our analysis in 
Sec. 2.1, we have used the average of these two sets of 6c values: our adopted 
values appear in the last column of the table. 

Two more recent calculations provide a valuable check that these Sc values 
are not suffering from severe systematic effects. Sagawa, van Giai and Suzuki ^-'^ 
have added RPA correlations to a Hartree-Fock calculation that incorporates 
charge-symmetry and charge-independence breaking forces in the mean-field 
potential to take account of isospin impurity in the core; the correlations, 



Table 4: Calculated isospin symmetry- breaking correction, 5c, in percent units. 



TH" OB'' SVS'^ NBO'' Adopted'^ 

Value 



0.18 0.15 0.00 0.12 0.17(3 

i^'O 0.32 0.15 0.29 0.24(3 

26'»A1 0.27 0.30 0.27 0.29(3 

3^01 0.64 0.57 0.33 0.61(3 

■'^^'"K 0.62 0.59 0.33 0.61(3 

42Sc 0.49 0.42 0.44 0.46(3 

^'^V 0.43 0.38 0.41(3 

5"Mn 0.51 0.35 0.43(3 

54Co 0.61 0.44 0.49 0.53(3 



Ref. * ; Ref. " ; " Ref. ''^ : Ref. 



" Average of OB and TH; assigned uncertainties reflect the 
relative scatter between these calculations. 



in essence, introduce a coupling to the isovector monopole giant resonance. 
The calculation is not constrained, however, to reproduce known separation 
energies. Finally, a large shell-model calculation has been mounted for the 
A = 10 case by Navratil, Barrett and Ormand^^. Both of these two new works 
have produced values of 6c very similar to, but actually smaller than those 
used in our analysis, i.e. worsening rather than helping the unitarity problem. 

The typical value of 5c is of order 0.4%. If the unitarity problem is to be 
solved by improvements in 5c, then all 5c values must be increased so that 
the typical value is raised to around 0.7%. There is no evidence whatsoever 
for such a shift from recent works. 



4 Speculative suggestions to resolve unitarity problem 

Since it is concluded in Sec. 3 that reasonable adjustments to the theoretical 
corrections, Sr, and Sc, are unlikely to resolve the unitarity problem posed 
by the nuclear result for Vud, we turn in this section to some more speculative 
alternatives. We discuss four suggestions, two that do not require extensions 
to the Standard Model and two that do. All but one require the introduction 
of important new physics to explain the apparent discrepancy. 



4-1 Saito-Thomas correction 



One suggestion to resolve the unitarity problem was made some years ago by 
Thomas It is based on a quark- meson coupling model, in which nuclear 
matter consists of non-overlapping nucleon (MIT) bags bound by the self- 
consistent exchange of a and w mesons in the mean-field approximation. The 
model is extended to include an isovcctor-vcctor meson (p) and an isovector- 
scalar meson [6). The coupled equations to be solved are very similar to those 
of Quantum Hadrodynamics but involve boundary conditions at the bag 
radius. As a consequence of a coupling to meson fields the quark mass in a 
medium becomes an effective one, 

m* =mi- {V„ ± \V5), i = u,d (24) 

with the upper sign for the up quark. Here and Vs are strengths of a- 
meson and (5-meson mean fields. At the quark level, the conserved vector 
current (CVC) hypothesis is broken if the up and down quarks have different 
masses, but for a free nucleon this breaking is second order in (m^ — m„). In a 
nuclear medium, however, this breaking may become first order in (m^ — m* ), 
but a critical constraint of the calculation is to show it reverts to (m^ — m„)^ 
at zero density. 

Saito-Thomas (ST) solve the coupled equations in an infinite nuclear- 
matter approximation, the solutions being functions of the matter density, 
Pb- Let ^i/j be the wavefunction of quark i bound in nucleon j in nuclear 
matter; then, the following overlap integral can be defined between quark i 
bound in a proton and quark i' bound in a neutron: 

UApb)^ I dV^^^i^e/n- (25) 

./Bag 

Calculations indicate that IwipB) varies linearly in ps, being unity at zero 
density. What is required for beta decay is the product of three overlap inte- 
grals 

|4d|'x|/„„|2x|/dd|2 = l-,5^"-'^ (26) 

and this product also is approximately linear in the density. The results of the 
calculations^^ are 

^^"^■■k = b X (^^^ , (27) 

where po is the saturation density for equilibrium nuclear matter, and b is 
in the range (0.15 — 0.20)% for bag radii in the range (0.6 — 1.0) fm. Thus 



jquark^ at densities of po/2, lies in the range 0.075% to 0.10% while, at po, it 
lies between 0.15% and 0.20%. If such a correction were to be applied to the 
results from nuclear superallowed transitions, the corresponding value for 
would then be increased by these amounts. 

However, in their original work, the authors themselves admit'*'' that their 
results are merely qualitative, having been derived from a model that deals only 
with nuclear matter. More recently, they have also recognized that a differ- 
ent model actually produces a considerably smaller effect. Furthermore, there 
remains the question of whether a simple addition of the Saito- Thomas correc- 
tion, S'^'^'^^, to the 6c already computed in a nucleons-only calculation is the 
correct approach at all, since it may lead to double counting. In the miclcons- 
only case, certain parameters are adjusted to reproduce Coulomb observables 
such as the b- and c-coefRcients of the isobaric mass multiplet equation, and 
proton and neutron separation energies. Thus, the parameters become effective 
ones, which in some unquantifiable way, actually contain quark effects. 

In summary, the Saito-Thomas correction continues to be unspecified for 
finite nuclei and now appears for several reasons to be a much smaller effect 
than was originally proposed. It seems clear that this effect, if it is applicable 
at all, would be insuflacient in itself to restore unitarity. 

4-2 Ad hoc parameterized fit 

Wilkinson has suggested that data like those displayed in Fig. 1 might be 
better fitted by a quadratic function in Z, representing a further correction 
of speculative origin. He considers the value of this fitted function at Z = 
to be the appropriate JFt- value to take forward in the determination of Vud- 
Alternatively, he has also considered the option of disregarding the structure- 
dependent corrections entirely and then fitting the uncorrected /f-values with 
a similar function quadratic in Z. There is no physical justification whatsoever 
for incorporating any Z-dependent corrections beyond those already accounted 
for in 6r and Sc- Furthermore, there is no statistically significant indication 
from the J-'t-vahic data in Fig. 1 that one is required. Consequently, until such 
time as new Tt-va,\ue data can demonstrate a clear residual .Z-dependence, 
this suggested explanation for the nuclear unitarity result must be regarded as 
unsatisfactory. 

4-3 Right-hand currents 

Because experimental evidence at low energies favors maximal parity violation 
in weak interactions, that condition has been built into the Standard Model. 



However, one class of possible extensions to the Standard Model would re- 
store parity symmetry at higher energies through the introduction of additional 
heavy charged gauge bosons that are predominantly right-handed in character. 
Under certain specific conditions, such models could explain the nuclear results 
described in Sec. 2.1. 

For example, an extension known as the manifest left-right symmetric 
model leads to a revised form^^ for Eq. (8) which includes a mixing angle (: 

If, under these conditions, we require that V^^ must satisfy unitarity, viz 

V^d = '^-Vl- (29) 

then, with the !Ft value taken from Eq. (7), we can derive a value for the 
mixing angle of C = 0.0015 ± 0.0007. Within the context of the manifest left- 
right symmetric model, this is a fully satisfactory result; not surprisingly, it is 
two standard deviations away from the Standard Model's pure V — A value. 



4-4 Scalar interaction 

If the Standard Model of pure V —A weak interactions is extended instead to in- 
clude scalar and tensor interactions, then the expression for the beta spectrum 
shape would include an additional term coming from the scalar interaction, 
which is inversely proportional to the electron energy. Since the J-t values 
include an integral over the spectrum shape, the presence of a non-zero scalar 
interaction would be reflected by a correction to the J^t values that is inversely 
proportional to the decay energy of the parent nucleus. To test this possibility, 
we fit the data in Fig. 1 by the function 



J't = k[l + hFl{W-^)\, (30) 

where k and hp are parameters determined from the fit. The latter, hp, is 
known as the Fierz interference term; a non-zero value for this term would 
signal the presence of a scalar interaction. In Eq. (30), 7^ = (1 — (aZ)^), a 
is the fine-structure constant, Z the atomic number of the daughter nucleus, 
and (yV~^) is the value of 1/W averaged over the electron spectrum, where W 
is the electron energy in electron rest-mass units. The value obtained for the 
Fierz interference term from the data in Fig. 1 is 



hp 



-0.0027 ±0.0029. 



(31) 



A negative sign is expected*^ for a positron emitter. The 90% confidence level 
upper limit is 

|6ir|< 0.0075. (32) 

To proceed to Vud, we need a value of J^t extrapolated to the Z = limit. To 

this end, wc fit the nine values of j{W~'^) by a second-order polynomial in Z 
and use this polynomial to provide the value of j{W~^) at Z = 0. With k and 
bp taken from the fit, Eq. (31), we obtain 

^i(O) = 3067.5 ±2.9 

\Vud\= 0.9748 ± 0.0006 [Nuclear + bF] 

^ 1/2. = 0.9983 ± 0.0016. [Nuclear + bF] (33) 

i 

This result is now in accord with unitarity but, of course, requires the presence 
of a non-zero scalar coupling to achieve that goal. 

Wc close this section by noting that all four suggested explanations we 
have presented for the non-unitarity result in Eq. (11) are entirely speculative 
and should be considered with caution. 

5 Future experimental prospects 

It is evident from the foregoing discussion that the current world data on 
s-re tantalizingly close to producing a definitive result on unitarity. The 
nuclear measurements have achieved the highest experimental precision but 
they arc now constrained by theoretical uncertainties. The neutron and pion 
measurements are, as yet, experimentally less precise than the nuclear ones, but 
they arc free from one of the more important sources of theoretical uncertainty, 
Sc- All three classes of measiirements are now being extended and improved at 
a number of laboratories, and there are good prospects for considerably reduced 
error bars within a few years. In combination with the re-visitation of the Ke3 
decay these results should settle the uncertainty over 5c and determine 
whether the deviation from unitarity, apparent from the nuclear result, is real 
or not. This is ample argument to justify considerable experimental activity. 
At the same time, it should not be forgotten that the full impact of these 
experimental advances will be diluted until there are theoretical improvements 
in the correction - the dominant source of theoretical uncertainty in all 
cases. Ultimately, to make significant improvements in the unitarity test, there 
will have to be advances in both theory and experiment. 



Of the three classes of experiment, that focusing on pion decay is currently 
farthest from the precision required for a meaningful unitarity test. We noted 
in Sec. 2.3 that a considerable improvement is anticipated in the foreseeable fu- 
ture, but the result will not reach a precision higher than that already achieved 
for neutron decay. There is considerable optimism, however, that the neutron 
measurements themselves can be improved as new experiments with ultra-cold 
neutrons come to fruition. The outcome can be expected seriously to challenge 
the nuclear experiments for experimental precision, but will take at least a few 
more years to do so. 

As to the nuclear experiments, the nine superallowed transitions whose 
ft values arc known to within a fraction of a percent have been the subject 
of intense scrutiny for at least the past three decades. All except have 
the special advantage that the superallowed branch from each is by far the 
dominant transition in its decay (> 99%). This means that the branching 
ratio for the superallowed transition can be determined to high precision from 
relatively imprecise measurements of the other weak transitions, which can 
simply be subtracted from 100%. Given the quantity of careful measurements 
already published, are there reasonable prospects for significant improvements 
in these decay measurements in the near future? Given the uncertainty in 
the theoretical corrections, which experiments can shed the most light on the 
efficacy of these corrections? 

If we begin by accepting that it is valuable for experiment to be at least 
a factor of two more precise than theory, then an examination of the world 
data shows that the Q-values for i°C, "O, ^S'^Al and ^'^V, the half-lives of 
^°C, ^^Cl and ^s^x, and the branching ratio for can all bear improve- 
ment. Such improvements will soon be feasible. The Q-values will reach the 
required level (and more) as mass measurements with new on-line Penning 
traps become possible; half-lives will likely yield to measurements with higher 
statistics as high-intensity beams of separated isotopes are developed for the 
new radioactive-beam facilities; and, finally, an improved branching-ratio mea- 
surement on has already been made with Gammasphere and simply awaits 
analysis . 

Qualitative improvements will also come as we increase the number of 
superallowed emitters accessible to precision studies. The greatest attention 
recently has been paid to the = emitters with A > 62, since these nuclei 
are expected to be produced at new radioactive-beam facilities, and their cal- 
culated structure-dependent corrections, Sc - Sns, (sec Eq.(5)) are predicted 
to be large^'**^'^. They could then provide a valuable test of the accuracy of 6c 
calculations. It is likely, though, that the required precision will not be attain- 
able for some time to come. First, with current Penning-trap technology, these 



nuclei are too short lived < lOOms) for their masses to be measured with 

the required precision. Second, the decays of these nuclei are of higher energy 
and each must involve numerous weak Gamow- Teller transitions that together 
will account for significant intensity in addition to the superallowed transition 

Branching-ratio measurements will thus be very demanding, particularly 
with the limited intensities likely to be available initially for these rather ex- 
otic nuclei. Finally, lifetime measurements will often be similarly constrained 
by statistics, although not in every case. The half-life of ^^Rb was recently 
obtained with fully adequate precision. 

More accessible in the short term will be the Tz = —1 superallowed emit- 
ters with 18 < A < 38. There is good reason to explore them. For example, 
the calculated value ^ of the total structure-dependent correction, 6c - Sns, for 

decay, though smaller than the values expected for the heavier nuclei, is 
actually 1.13% about a factor of two larger than for any other case currently 
known - while ^^Mg has a rather low value of 0.51%. If such large differences 
are confirmed by the measured /f-values, then it will do much to increase 
our confidence in the calculated isospin symmetry-breaking corrections. To be 
sure, these decays will provide a challenge, particularly in the measurement 
of their branching ratios, but the required precision should be achievable with 
isotope-separated beams that are currently available. In fact, such experiments 
are already underway at the Texas A&M cyclotron. Results for the decay of 
^^Mg are currently being prepared for publication. 

6 Conclusions 

The current world data on superallowed 0+ — > 0+ beta decays lead to a 
self-consistent set of ^i-values that agree with CVC but differ provocatively, 
though not yet definitively, from the expectation of CKM unitarity. There are 
no evident defects in the calculated radiative and isospin symmetry-breaking 
corrections that could remove the problem, so, if any progress is to be made in 
firmly establishing (or eliminating) the discrepancy with unitarity, additional 
experiments are required. We have indicated what some relevant nuclear ex- 
periments might be. 

In the past decade there have been significant improvements in the mea- 
surements of the neutron lifetime and beta asymmetry, and further improve- 
ments are promised in the near future. It is likely that these studies in neutron 
decay will soon approach the results from nuclear superallowed decays in pre- 
cision; and when there, they will have the advantage in that their results are 
not dependent on nuclear-structure corrections. So far, though, the neutron 
result is quite consistent with the nuclear one. 



Clearly, there is strong motivation to pursue experiments, not only on the 

neutron and nuclear front but also in re- visiting the K(,3 decay, since, if firmly 
established, a discrepancy with unitarity would indicate the need for important 
new physics. 
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